Creation of new germplasm resources, development of SSR markers, and screening of monoterpene synthases in thyme

Background Thyme derived essential oil and its components have numerous applications in pharmaceutical, food, and cosmetic industries, owing to their antibacterial, antifungal, and antiviral properties. To obtain thyme essential oil with different terpene composition, we developed new germplasm resources using the conventional hybridization approach. Results Phenotypic characteristics, including essential oil yield and composition, glandular trichome density, plant type, and fertility, of three wild Chinese and seven European thyme species were evaluated. Male-sterile and male-fertile thyme species were crossed in different combinations, and two F1 populations derived from Thymus longicaulis (Tl) × T. vulgaris ‘Fragrantissimus’ (Tvf) and T. vulgaris ‘Elsbeth’ (Tve) × T. quinquecostatus (Tq) crosses were selected, with essential oil yield and terpene content as the main breeding goals. Simultaneously, simple sequence repeat (SSR) primers were developed based on the whole-genome sequence of T. quinquecostatus to authenticate the F1 hybrids. A total of 300 primer pairs were selected, and polymerase chain reaction (PCR) was carried out on the parents of the two hybrid populations (Tl, Tvf, Tve, and Tq). Based on the chemotype of the parents and their F1 progenies, we examined the expression of genes encoding two γ-terpinene synthases, one α-terpineol synthase, and maybe one geraniol synthase in all genotypes by quantitative real-time PCR (qRT-PCR). Conclusion We used hybridization to create new germplasm resources of thyme, developed SSR markers based on the whole-genome sequence of T. quinquecostatus, and screened the expression of monoterpene synthase genes in thyme. The results of this study provide a strong foundation for the creation of new germplasm resources, construction of the genetic linkage maps, and identification of quantitative trait loci (QTLs), and help gain insight into the mechanism of monoterpenoids biosynthesis in thyme. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-022-04029-2.

To date, a few studies have been conducted by scientific research institutions and companies on thyme breeding, especially in European countries, and few of the high-performance industrial and horticultural varieties developed in these studies have been cultivated. In the 4th International Symposium on the Breeding of Medicinal Aromatic Plants, a breeding study was reported, in which male-sterile and male-fertile clones were crossed to optimize the terpenoid content and yield of thyme (T. vulgaris) [12]. In this study, 56 new hybrids, which were tested by assessing the homogeneity, dry weight, essential oil yield, winter frost tolerance, and seed production potential of the parents, were obtained from 2000 to 2002. The most dominant hybrid, named 'Varico 3' , showed 4.9% essential oil yield and thymol-type chemotype [12]. Such hybrid varieties have been entered into the market. Clone T-12, with high phenol content, was selected from among 10 T. vulgaris clones [13]. The global collection of Thymus resources is highly diverse, with more than 300 species, which are native to the Mediterranean basin and are widely distributed in the temperate regions of Europe, North Africa, and Asia [14]. In China, Thymus species are mainly distributed in the northwest, north, and northeast regions based on the Flora of China [15]. Therefore, according to the breeding goals, wild Chinese thyme species could be crossed with European thyme species to develop a series of new varieties for applications in different fields.
With the rapid development of modern molecular biology approaches and genome sequencing technologies, DNA-based molecular markers have become an important tool for cultivar identification, fingerprinting [16,17], and genetic diversity analysis [18][19][20]. To improve the industrial applications of medicinal and aromatic plants, breeders often select plants with high genetic divergence and essential oil content [17]. Phenotypic variation could be very valuable for molecular breeding approaches such as marker-assisted selection (MAS), which has been very helpful in elucidating the genetic diversity of plant species. A dendrogram, based on cluster analysis, showed that T. daenensis and T. fallax are clearly distinct from the other Thymus species, indicating that T. daenensis shares some genetic similarity with T. fallax [21]. Previously, several studies used randomly amplified polymorphic DNA (RAPD) markers to investigate the genetic diversity and essential oil composition of various Thymus species as well as the phylogenetic relationship among these species [11,17,22].
In this study, male-sterile and male-fertile thyme varieties were crossed in two different combinations, T. longicaulis (Tl) × T. vulgaris 'Fragrantissimus' (Tvf ) and T. vulgaris 'Elsbeth' (Tve) × T. quinquecostatus (Tq), to generate two F 1 populations. Then, simple sequence repeat (SSR) markers were developed based on the whole-genome sequence of T. quinquecostatus [36], to authenticate all F 1 individuals. In addition, the expression profiles of two γ-terpinene synthase genes (Tq13G005250.1 and Tq02G002290.1), one geraniol synthase gene (Tq04G005190.1), and maybe one α-terpineol synthase gene (Tq03G001560.1) were analyzed by quantitative real-time PCR (qRT-PCR). Overall, this study provides a valuable collection of new thyme varieties, which could be used for MAS and the verification of TPS gene function in future studies.

Phenotypic evaluation of different thyme species and construction of F 1 hybrid populations
Thyme is a herbaceous perennial or sub-shrub with valuable medicinal and aromatic properties. Thyme plants are of two types, depending on their growth habit: erecttype and creeping-type. These two plant types display remarkable differences in morphology. The plant types of 10 different thyme species are shown in Fig. 1a Fig. 1b. Tr, Tve, Tg, Tt, Ts, and Tl contained only stigma and no pollen (stamen), and were therefore male-sterile. By contrast, Tq, Tqp, Tm, and Tvf possessed both stigma and pollen (stamen), and therefore were categorized as male-fertile ( Table 1). The increased genetic diversity of thyme species could be attributed partially to the consistent introgression of wild Chinese thyme germplasm into the male-sterile and erect-type European thyme germplasm during long cultivation periods, and partially to the adaptation of thyme species to new environments in new geographical locations.
Glandular trichomes are specialized hairs that originate from the epidermal cells of flowers, leaves, and stems. These organs exist as two types, peltate and capitate, on the surface of approximately 30% of all vascular plants, including lavender, thyme, rosemary, oregano, basil, and other Lamiaceae species [37]. Glandular trichomes are responsible for a significant portion of the secondary metabolite of a plant [38], and serve as the storage and synthesis sites of terpenoids [39]. The regulation of glandular trichome formation related genes potentially underlies the regulation of glandular trichome density (number of glandular trichomes per unit area) for increasing the terpenoid content of plants.
Tve and Tr showed the highest glandular trichome density on both the adaxial and abaxial leaf surfaces (Fig. 1c, e). Glandular trichome density on the adaxial leaf surface was the second highest in Tvf and Tt. The glandular trichome density of leaves, overall, decreased in the following order: Tve, Tr, Tvf, Tt, and Ts. The essential oil yield of thyme species is shown in Fig. 1d, f and Supplementary   Fig. S1. Tvf showed the highest essential oil yield (1.25 mL·100 g − 1 ), followed by Tt (1.20 mL·100 g − 1 ), Tr (0.95 mL·100 g − 1 ), Tve (0.75 mL·100 g − 1 ), and Tl (0.70 mL·100 g − 1 ). The essential oil yield of Tq, Tqp, Tm, Tg, and Ts varied between 0.40 mL·100 g − 1 and 0.50 mL·100 g − 1 ( Table 1; Supplementary Fig. S1). There was a certain correlation between essential oil yield and glandular trichome density, the higher the density of glandular trichome, the higher the yield of essential oil (Fig. 1e, f ).
Cluster analysis of the 20 main compounds found in the essential oil of 10 different thyme species (Fig. 2b) revealed four clusters. Tq and Tqp, which contained carvacrol as the most abundant compound (carvacroltype essential oil), with relative contents of 20.74% and 48.37%, respectively, clustered together; Tm, Tve, Tr, Tt, Tg, and Ts grouped together, and contained thymol as the most abundant compound (thymol-type essential oil); and Tl and Tvf clustered separately (geraniol-type and α-terpineol-type essential oil, respectively). Principal component analysis (PCA) analysis was carried out on the main compounds found in the essential oils of all 10 thyme species (Fig. 2c, d). The results showed that the relatively high contents of thymol, carvacrol, p-cymene, and γ-terpinene contributed greatly to the volatile components of the essential oil of all 10 species. Tq and Tqp grouped together in the first quadrant, and their corresponding characteristic volatile substances included carvacrol, p-cymene, and γ-terpinene; Tvf and Tl were distributed in the third quadrant, which corresponded to volatiles geraniol and α-terpineol; Tm, Tve, Tr, Tt, Tg, and Ts were distributed in the fourth quadrant, and the corresponding predominant volatile was thymol (Fig. 2d); these results verified the results of cluster analysis.

Development and application of SSR markers
The chromosome-level genome assembly and annotation using high-fidelity (HiFi) and chromatin conformation capture (Hi-C) technologies revealed 13 chromosomes in T. quinquecostatus, with a total length of 528.66 Mb, 70.61% (373.28 Mb) of which was annotated as highly repetitive [36].  S2c). Among the dinucleotide repeats, the CT/AG-type repeat was the most abundant (23.   AAAT/ATTT (0.40%) was the only tetranucleotide repeat type identified in the T. quinquecostatus genome ( Supplementary Fig. S2c). The length of SSR loci ranged from 18 to 87 bp in the T. quinquecostatus genome, with 10 bp SSRs being the most abundant (61,865, accounting for 32.00% of all SSR loci) and 25 bp SSR loci being the least abundant (1,356, 0.71%). The number of SSRs gradually decreased with the increase in repeat length (Supplementary Fig. S2d).
Male-sterile (without pollen) and male-fertile (with pollen) thyme species were crossed as female and male parents, respectively, in different combinations. Finally, two F 1 hybrid populations were obtained from two crosses: T. longicaulis × T. vulgaris 'Fragrantissimus' (Tl × Tvf, 14 lines) and T. vulgaris 'Elsbeth' × T. quinquecostatus (Tve × Tq, 11 lines) (Supplementary Table S1). To design SSR markers for the verification of F 1 progenies, 300 primer pairs were screened by performing PCR amplification on the parental lines of the two crosses (Supplementary Table  S2). Analysis of the PCR products revealed 1-2 polymorphic bands between the two parents of each cross. After many repetitions, the primers showing clear and stable banding patterns were selected (Fig. 3a, d).
Eighteen SSR markers were co-dominant in the Tl × Tvf population (Supplementary Table S3). For example, TqSSR289 amplified band 1 in the female parent T1 and band 2 in the male parent Tvf (Fig. 3a); TqSSR292 amplified band 3 in the female parent T1 and band 4 in the male parent Tvf (Fig. 3a). Similarly, 23 SSR markers were co-dominant in the Tve × Tq population (Supplementary Table S4). For example, TqSSR284 amplified band 1 in the female parent Tve and band 2 in the male parent Tq (Fig. 3d), so as to be used for the identification and verification of the hybrid progenies of the combination, accounting for 7.60% of 300 primer pairs. These co-dominant SSR primers were used to identify F 1 individuals in the two hybrid populations. Progenies with complementary parental bands or only paternal-specific bands were true hybrids, and those with only maternal-specific bands were pseudo-hybrids or inbreds. Based on the genotyping results, 14 lines in the Tl × Tvf progeny (Fig. 3b, c), and 11 lines in the Tve × Tq progeny were identified as true hybrids (Fig. 3e, f ).

Phenotypic evaluation of thyme species and construction of two F 1 hybrid populations
In this study, the essential oil yield of erect-type thyme was higher than that of creeping-type thyme, and there was a positive correlation between the density of glandular trichomes and essential oil yield (Fig. 1e, f ). Some wild Chinese thyme species showed a long flowering period of up to four months, while some European thyme species showed a short flowering period of one month ( Table 1). The tillering capacity of creeping-type species was greater than that of erect-type species. In addition, the harvest time, extraction method, and stem/leaf ratio of thyme species likely affected their essential oil yield. In this study, the type and relative content of VOCs in the essential oil varied significantly among the 10 different thyme species, which further illustrated the diversity of the Thymus genus (Table 2). A previous study reported that the essential oil composition of a given thyme species and the contents of individual VOCs vary among different environments [40]. Additionally, extraction methods affect the composition and content of essential oil [17]. Therefore, the effects of species, region, and extraction method should be considered when studying the content and composition of thyme essential oil (Supplementary Fig. S1; Table 2). A certain correlation was also detected between the color and composition of thyme essential oil. Essential oil extracted from α-terpineol-type thyme was milky white whereas that extracted from geranioltype thyme was light yellow. Essential oils extracted from thymol-type and carvacrol-type thyme species were golden yellow. The most well-studied volatile compounds of Chinese and European thyme species are thymol and carvacrol, since these two compounds are the two most abundant VOCs in most thyme species [4,8,9]. In this study, the chemical types of the seven European thyme species were mainly thymol-, geraniol-, and α-terpineoltype, while those of the three wild Chinese thyme species were carvacrol-and thymol-type. The phenomenon of male sterility in plants enables cross breeding while avoiding complicated process of emasculation. Darwin reported that thyme plants in southern England produce two kinds of flowers: one with intact male and female organs, and the other (a small flower) with no or completely cut anthers; the latter are completely male-sterile [41]. In different ecological environments, the proportion of female plants in different thyme species is more than 50% [42]. Because thyme flowers are very small in size, emasculation and crossing are difficult. Therefore, if some thyme plants possess female reproductive organs, crossing is very easy. Among the 10 thyme species investigated in this study, Tr, Tve, Tg, Tt, Ts, and Tl were male-sterile, while Tq, Tqp, Tm, and Tvf were male-fertile. Therefore, Tve and T1 were used as female parents, and Tq and Tvf were used as male parents to generate two hybrid thyme populations. Male sterility is jointly determined by genetic and environmental factors. The male-sterile thyme lines used in this study were introduced from Europe; therefore, the phenomenon of male sterility in these lines may be caused by maladjustment to the new environment.
Chinese thyme plants are widely distributed in the northwest, north, and northeast regions in China [15]. Therefore, the research on Chinese thyme species should be intensified, and wild Chinese thyme species should be crossed with European thyme species to breed new species suitable for ornamental, medicinal, and edible applications to promote the development of the thyme planting industry. Hybrid breeding is an important means of germplasm resource innovation [12]. The advantages of both parents can be introgressed into a single genetic background through interspecific hybridization. In this study, new thyme germplasm resources were generated by crossing different chemical types and plant types of thyme.

SSR marker development and application
Hybrid breeding is an important means of germplasm resource innovation. The advantages of parents can be combined in a single genetic background through interspecific hybridization. In this study, new thyme germplasm resources were generated by hybridizing different chemical types and plant types. SSR markers are simple, time-saving, cost-effective, reproducible, and stable, and therefore have been widely used for the identification of the hybrid progenies in different species [16,43]. Early identification and selection of hybrid progeny is an important link in cross breeding. Therefore, molecular marker technology, combined with morphological observation, can be used to identify hybrid progenies with high accuracy.
The SSR-based genotyping process first involves the screening of co-dominant SSR primers in parents, followed by the identification of progenies [43]. However, there may be errors in SSR marker-based identification of hybrid progenies using a single primer pair. For example, in this study, SSR primers TqSSR119 and TqSSR124 were used to identify F 1 progenies derived from the Tl × Tvf cross. Analysis of the same sample using different primer pairs can produce different results, possibly because of the high heterozygosity of thyme, lead to some of the parents can't dominance or appear some non-parental bands. This may also be caused by the process of meiotic division during gamete formation; exchange of DNA between homologous chromosomes during recombination at marker loci can lead to the disappearance of bands in progenies.
Generally, the phenotype of F 1 progenies is intermediate between that of the two parents or biased towards the phenotype of the female or male parent, which also provides the possibility of breeding excellent new varieties [43]. In this study, among the 14 Tl × Tvf F 1 lines, eight lines showed the same VOC profiles in leaves as their female parent Tl; one line was similar to the male parent Tvf; and five lines were dissimilar to both parents. This conclusion was supported by the 11 F 1 lines derived from the Tve × Tq cross; five of these F 1 lines showed the same leaf VOC profiles as their female parent Tve, and the remaining six lines were different from both their parents. However, some F 1 progenies also showed some volatile compounds that were superior to the both parents or not found in either parent; for example, thymol and carvacrol were found in the progeny of the Tl × Tvf cross, and α-terpineol and α-terpineol acetate were found in the Tve × Tq population. These results suggest that new varieties with related VOC profiles could be developed in thyme through cross breeding. Heterosis may also be present, and the greater the phenotypic differences between parents, the stronger is the heterosis expected to be in their progeny.

Bioinformatics analysis and screening of TPSs in thyme
Monoterpenoid biosynthesis begins with GPP, which is the precursor of all monoterpenoids, and yields α-terpinyl cations, which are highly unstable intermediates that can then be converted to specific monoterpenoids by certain monoterpene synthases such as γ-terpinene and α-terpineol [44]. GPP also serves as the synthetic precursor of geraniol, linalool, myrcene, and ocimene, which are formed through catalysis by different monoterpene synthases. In addition, cytochrome P450 monooxygenase 71D (CYP71D) proteins subfamily and short-chain dehydrogenases/reductases (SDRs) are involved in further modification of the γ-terpinene framework to produce p-cymene, thymol, and carvacrol [36]. The synthetic pathway of some monoterpenes is shown in Fig. 6b. Functions of the γ-terpinene synthase gene TcTPS02 and α-terpineol synthase gene TcTPS05 were previously validated in T. caespititius [29]. Similarly, functions of the γ-terpinene synthase gene TvTPS2 and three cytochrome P450 genes (TvCYP71D179, TvCYP71D180, and TvCYP71D507) were validated in T. vulgaris [45,46].
The results of qRT-PCR analysis showed that the screened TPSs were expressed in hybrid thyme progenies, and these results were consistent with the relative contents of the TPS-catalyzed products in thyme. The expression levels of two γ-terpinene synthase genes Tq02G002290.1 and Tq13G005250.1 were verified in the progeny derived from the cross between Tve and Tq, whose chemotypes were thymol-and carvacroltype, respectively. Similarly, the expression levels of α-terpineol synthase gene Tq03G001560.1 and geraniol synthase gene Tq04G005190.1 were validated in the progeny of Tl and Tvf, whose chemotypes were α-terpineoland geraniol-type, respectively. Additionally, the relative expression levels of Tq02G002290.1, Tq13G005250.1, Tq03G001560.1, and Tq04G005190.1 were consistent with the relative contents of catalytic products in some F 1 lines. Therefore, this method enables only a preliminary screening of gene function and provides a basis for further gene function verification. The above results lay a strong foundation for the creation of new germplasm resources, construction of the genetic linkage maps, mapping of quantitative trait loci (QTLs), and MAS, and provide insight into the mechanism of monoterpenoids biosynthesis in thyme.

Conclusion
Thyme is a multi-purpose plant with a wide range of applications in the pharmaceutical, food, and cosmetic industries, owing to its high essential oil content. To obtain thyme essential oil with different terpene compositions, we cross wild Chinese thyme species with European thyme species, and new germplasm resources have developed using the conventional hybridization approach. Two F 1 populations were obtained, simultaneously, SSR primers were developed based on the whole-genome sequence of T. quinquecostatus to authenticate the F 1 hybrids. Based on the chemotype of the parents and their F 1 progenies, we examined the expression of genes encoding two γ-terpinene synthases, one α-terpineol synthase, and maybe one geraniol synthase in all genotypes by qRT-PCR. The results of this study provide a strong foundation for the creation of new germplasm resources, construction of the genetic linkage maps, and identification of QTLs responsible for the terpene compositions and essential oil. We screened monoterpene synthases in thyme to gain insight into the mechanism of monoterpenoids biosynthesis in thyme.

Hybrid breeding design
Male-sterile (without pollen) thyme species were crossed as the female parent with male-fertile (with pollen) thyme species, with the essential oil composition and yield as the main breeding goals. Different cross combinations were designed, and the F 1 hybrid combinations of thyme were constructed by cross breeding in 2020. Finally, two F 1 populations derived from T. longicaulis (Tl) × T. vulgaris 'Fragrantissimus' (Tvf ) and T. vulgaris 'Elsbeth' (Tve) × T. quinquecostatus (Tq) crosses were selected for further analysis.

DNA extraction
Leaves were collected from the parents and F 1 progenies of Tl × Tvf and Tve × Tq crosses, immediately frozen in liquid nitrogen, and stored at -80 °C. DNA was extracted from the frozen leaf samples using the DNA Secure Plant Kit (Tiangen, China). DNA concentration and quality were assessed by 1% agarose gel electrophoresis and with a 2.0 Fluorometer (Life Technologies, CA, USA).

RNA extraction and cDNA synthesis
Leaves of T. longicaulis, T. vulgaris 'Fragrantissimus' , T. vulgaris 'Elsbeth' , T. quinquecostatus, and their F 1 lines were frozen in liquid nitrogen and stored at -80 ℃. Total RNA was extracted from the frozen leaves using the Easy Spin RNA extraction kit (Sangon Biotech, Shanghai, China). The isolated total RNA was treated with DNase I, and then purified with the RNA clean kit (Promega, Madison, WI, USA). The concentration of each RNA sample was determined using NanoDrop spectrophotometer (Thermo Fisher Scientific Inc., USA) and 2.0 Fluorometer (Life Technologies, CA, USA). RNA integrity was analyzed using Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA). Then, cDNA was synthesized using the HiScript Reverse Transcriptase Kit (Vazyme, China), according to the manufacturer's instructions [33].

Essential oil extraction
The essential oil of 10 thyme species was isolated by steam distillation at 180-200 ℃ for 90 min. Essential oil yield (%) was calculated as volume (ml) of the isolated oil per 100 g of dry plant material. The isolated essential oil was dried using anhydrous sodium sulfate, and stored at 4 ℃ until needed for further analysis [11].

Analysis of essential oil composition
The essential oil composition of 10 thyme species was analyzed by gas chromatography-mass spectrometry (GC-MS) using Agilent 7890 A-7000B gas chromatograph (Agilent, USA), equipped with Agilent 5975 C MS detector (Agilent, USA). Using the HP-5MS (30 m, 250 m ID, 0.25 μm film thickness) capillary column, volatiles were separated using the following temperature program: 5 min at 60 °C; increased to 220 °C at the rate of 4 °C/min; increased to 250 °C at the rate of 60 °C/min; hold at 250 °C for 5 min. The following parameters were used: injector and detector temperature, 250 °C; carrier gas, He; flow rate, 1 m/min; split ratio, 1:10; acquisition range, 50-500 m/z in electronimpact mode; ionization voltage, 70 eV; and injected sample volume, 1 µl. The determination of the content of each compound (%) was based on the normalization of GC peak areas. The identification of essential oil components was based on the comparison of retention indices (RIs), relative to a homologous series of n-alkanes (C7-C40), and mass spectra (MS) from the NIST (v14.0) library and data from scientific literature [48]. RIs were based on the equation: where RT(x), RT(z), and RT(z + 1) for the composition, and the number of carbons Z and Z + 1 for the retention time of the normal alkane.

Analysis of leaf VOC profiles of thyme
The leaf VOC profiles of T. vulgaris 'Elsbeth' , T. quinquecostatus, T. longicaulis, T. vulgaris 'Fragrantissimus' , and their F 1 lines were detected via headspace solid-phase microextraction (HS-SPME). Briefly, 0. sealed using crimp-top caps with TFE-silicone headspace septa (Agilent, Palo Alto, CA, USA). Subsequently, each vial was immediately incubated at 40℃ for 30 min. Then, to absorb the volatiles, the headspace of each vial was exposed to 100 μm coating fiber polydimethylsiloxane (Supelco, Inc., Bellefonte, PA, USA) for 30 min. All VOCs on the coating fiber were analyzed by GC-MS using Model 7890 A GC instrument and 7000B mass spectrometer (Agilent, Palo Alto, CA, USA) [49].
The GC-MS conditions were as follows: injector temperature, 250 ℃; transfer line temperature, 250 ℃, respectively; column temperature, initially maintained at 50 ℃ for 3 min, gradually increased to 150 ℃ at 4 ℃/ min for 2 min, and finally raised to 250 ℃ at 8 ℃/min for 5 min; carrier gas (helium) flow rate, 1 ml/min; injection, splitless mode; ionization voltage, 70 eV; source temperature, 250 ℃; and MS range, 35-500 m/z. Agilent Mass-Hunter 5.0 was used to analyze the chromatograms and MS. VOCs were identified by comparing the retention times of individual peaks with those of authentic standards, and MS were determined based on the NIST v14.0 MS database and the data reported previously [50]. RIs were calculated using the following equation: where RT(x), RT(z) and RT(z + 1) for the composition, and the number of carbons Z and Z + 1 for the retention time of the normal alkane.

Density of glandular trichomes
Glandular trichomes were visualized using a stereomicroscope (Leica DVM6, Germany). The number of glandular trichomes within a certain leaf area was counted using the ImageJ software. The average glandular trichome density was calculated based on three plants.

Statistical analysis
Data were expressed as the mean ± standard deviation of three biological replicates. Statistical analysis, including variance analysis, hierarchical clustering analysis, and correlation analysis, was performed using IBM SPSS Statistics for Windows, version 19.0 (Armonk, USA). Significant differences among the different genotypes were tested using a one-way analysis of variance (ANOVA), followed by Duncan's multiple range test at 5% probability level (p ≤ 0.05).

DMAPP
Dimethylallyl diphosphate EO Essential oil FPP Farnesyl diphosphate GC Gas chromatography